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ABSTRACT 

Cool stars at giant and supergiant evolutionary phases present low velocity and high 
density winds, responsible for the observed high mass-loss rates. Although presenting 
high luminosities, radiation pressure on dust particles is not sufficient to explain the 
wind acceleration process. Among the possible solutions to this still unsolved problem, 
Alfven waves are, probably, the most interesting for their high efficiency in transfering 
energy and momentum to the wind. Typically, models of Alfven wave driven winds 
result in high velocity winds if they are not highly damped. In this work we determine 
self-consistently the magnetic field geometry and solve the momentum, energy and 
mass conservation equations, to demonstrate that even a low damped Alfven wave 
flux is able to reproduce the low velocity wind. We show that the magnetic fluxtubes 
expand with a super-radial factor S > 30 near the stellar surface, larger than that 
used in previous semi-empirical models. The rapid expansion results in a strong spa- 
tial dilution of the wave flux. We obtained the wind parameter profiles for a typical 
supergiant star of 16 M@, The wind is accelerated in a narrow region, coincident with 
the region of high divergence of the magnetic field lines, up to 100 km s _1 . For the 
temperature, we obtained a slight decrease near the surface for low damped waves, 
because the wave heating mechanism is less effective than the radiative losses. The 
peak temperature occurs at r ~ 1.5 ro reaching 6000 K. Propagating outwards, the 
wind cools down mainly due to adiabatic expansion. 
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1 INTRODUCTION 

Cool giant and supergiant stars are known to present con- 
tinuous mass loss process occuring at high rates, typically 
1CT 10 - 10~ 5 M yr"\ but in low velocity winds (u™ < 300 
km s" 1 ) (Dupree 1986, Lamers & Cassinelli 1999). After 
decades of theoretical and observational studies, the mecha- 
nisms in which the wind acceleration occurs are still poorly 
understood. Due to the stellar high luminosity and low effec- 
tive temperature several authors have proposed radiativelly 
dust driven models to explain the observed wind proper- 
ties (Liberatore, Lafon & Berruyer 2001, Elitzur & Ivezic 
2001, Woitke & Niccolini 2005). At the pulsating phase, 
acoustic waves generate high density shells that could al- 
low dust to form near the stellar surface. In this scenario, 
radiative pressure on grains transfers momentum to these 
particles being responsible for their acceleration and, if gas 
and dust are dynamically well coupled, grains drag the gas 
outwards resulting in the mass ejection. However, at station- 
ary envelopes (e.g. pre-AGB phase) dust driven theoretical 
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models have lately failed in reproducing the wind proper- 
ties, mainly because the dust-gas coupling is not effective 
(Sandin & Hofner 2003). Observationally, Guandalini et al. 
(2005) found no strong correlation between the mass loss 
rates and the luminosities of AGB stars. Their main conclu- 
sion is that, if radiative pressure is important in powering 
these stellar winds, it must occurs in addition to other mech- 
anism. Another drawback to the radiation pressure models 
is the need for the dust formation region to be close to the 
star. Recent high resolution Doppler measurements show 
that winds are mainly accelerated near the stellar surface 
(r < 1.3 R*) (Airapetian, Carpenter & Ofman 2003), while 
grains are expected to grow and survive at even larger dis- 
tances. 

In this sense, another mechanism must be used to accel- 
erate the gas near surface. The most promising mechanism 
for the winds of cool stars is the transfer of momentum and 
energy to the wind from MHD waves. Alfven waves driven 
wind models are known to result in high velocity winds, 
unlike what is measured for giant cool stars, because these 
MHD waves are, in general, weakly damped. Hartmann & 
MacGregor (1980) showed that it would be possible to re- 
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produce the observed low wind velocities and the high mass 
loss rates of the cool giant and supergiant stars if some kind 
of wave damping mechanism is effective at the wind basis 
(r < 2 R»). In that work, they confirmed this assumption 
using the ion-friction damping, which has a damping length 
proportional to P 2 , where P is the wave period. Holzer, Fla 
& Leer (1983) questioned this result in terms of an unnatural 
fine-tuning for the wave flux period, and argued that stars 
would, more reasonably, present a wide variety of wave fre- 
quencies depending on the generation mechanisms, and for 
most of them the damping would be uneffective. However, 
Jatenco-Pereira & Opher (1989) studied the effects of both 
different damping mechanisms and the magnetic field diver- 
gence, showing that the fine-tuning of the wave period is un- 
necessary since there are several other damping mechanisms 
that would act on the frequency spectrum. Also, they showed 
that the magnetic field divergent geometry can rapidly di- 
lute the wave flux and also slow down the wind. Their mag- 
netic field geometry was based on empirical relations found 
from observations of the solar wind. This because in gen- 
eral, in lack of direct measurements of the magnetic field 
fluctuations and structure in other stars, we have to simply 
extrapolate our knowledge from solar observations. 

Parker (1958) proposed that the solar wind is acceler- 
ated by the strong thermal pressure gradient in the transi- 
tion region, where the gas temperature increases from ~ 10 4 
K up to 10 6 K. However, it became clear in the following 
years that neither the solar radiation, nor acoustic waves 
generated in the photosphere, could account for the heat- 
ing of the coronal base. The origin of the energy respon- 
sible for the heating and acceleration of the plasma is be- 
lieved to be from both magnetic field reconnections above 
the photosphere (Axford & McKenzie 1996) and convective 
motions under the stellar surface. The convective energy is 
transferred up to the atmosphere from the perturbations 
generated in the field line footpoints (Cranmer & van Bal- 
legooijen 2005). These perturbations propagate outwards as 
Alfven waves, heating the gas as they are damped (Suzuki 
& Inutsuka 2005) . Regarding the magnetic field lines, Holzer 
& Leer (1980) and Jatenco-Pereira & Opher (1989) realized 
that the super-radial geometry of the magnetic field funnels 
at the coronal holes could have a significant impact on the 
mass flux and wind speed. Esser et al. (2005), showed that 
the models considering highly diverging magnetic funnels 
explain better the observed data for the Sun. Also, Tu et 
al. (2005) established that the solar wind acceleration initi- 
ates in the magnetic funnels at heights lower than 2.10 9 cm, 
coincident with the high divergence region. 

In this work, wc model the acceleration and heating of 
a late-type supergiant stellar wind considering an outward 
flux of Alfven waves. We solve the MHD equations to, self- 
consistently, determine the magnetic field geometry and the 
wind temperature, density and velocity profiles. In section 
2, we describe the model basic equations. In section 3 we 
present the results and the discussions, followed, by the work 
conclusions. 



2 THE MODEL 

The wind equations are based on mass, momentum, energy 
and magnetic flux conservation. The first is given by: 



puA(r) = p u A(r ) 



(1) 



where u is the flow velocity, p is the gas density and A(r) 
is the flow cross-section area at a distance r from the center 
of the star. The index "0" indicates the variable is being 
evaluated at the stellar surface (r = r ). 

Assuming a steady flow, the momentum equation can 
be written as: 
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where P = pk B T/m is the thermal pressure, k B is the 
Boltzmann constant, T is the gas temperature, m is the 
mean mass per particle, G the gravitational constant and 
SB the wave magnetic field amplitude. In Equation (2), the 
right hand side contains the gravitational force density and 
the thermal and wave pressure gradients, respectively. The 
last two terms represent the Lorentz force. The wave am- 
plitude (SB) is related to the wave energy density (e) by 
e=((SBf)/(^). 



2.1 Thin fluxtube approximation 

Typically, considering magnetic field strengths > 1 G, the 
wind basis is characterized by the relation B 2 /8tt 3> P > 
pu 2 /2, i.e. the plasma is magnetically dominated. In this 
case, if we assume the wind to be initiated at funnels an- 
chored at the stellar surface, which are surrounded by a 
plasma with lower magnetic field strength, the magnetic 
pressure inside will push the gas and the funnel field lines 
will expand. The funnel cross-section radius (TV) will grow 
super-radially up to a limit value (1Z m ). This limiting cross- 
section radius depends both on the relation between external 
and internal magnetic field strengths and on the filling fac- 
tor (q). As the area increases, the internal magnetic strength 
diminishes until the equilibrium between internal and exter- 
nal magnetic pressures is reached. If the internal magnetic 
field strength is much larger than the external, the flux tubes 
cross-section will depend on the filling factor only. The fill- 
ing factor is the ratio between the area of the stellar surface 
covered by funnels and the total area. The averaged maxi- 
mum area that a funnel could reach would be A m — A(ro) /a 
or, in terms of the cross section radius: 
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For the quiet Sun, the funnels that merge to form the coronal 
holes cover about 10% of the total surface. 

To evaluate the tube expansion at the wind basis, we as- 
sume the plasma to be magnetically dominated and the left 
hand side of Equation (2) may be neglected if compared to 
the other terms. Then, by using V • B = and a power series 
expansion method proposed by Pneuman, Solanki & Stenflo 
(1986), we can determine self-consistently the magnetic field 
geometry without assuming any empirical function for the 
funnel cross-section with distance. 

Following Pneuman, Solanki & Stenflo (1986), using the 
thin fluxtube aproximation, Equation (2) and V • B — arc 
described near stellar surface by: 
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where H = ksTr 2 /GmM* is the scale height. 

Expanding all variables as power series in y (i.e. along 
the tube radius) , and neglecting terms of orders higher than 
2, Equations (4) — (6) give rise to a differential equation for 
the fluxtube cross section: 



(7) 
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where (3 = 47rP(r )/ B 2 (r ), ql is the filling factor and B ex t 
is the magnetic field strength outside the fluxtube. In the 
following calculations we fixed its value to be 10~ 3 £>(r). 

To simplify the set of equations, we will define the funnel 
area expansion as a function of radial distance by: 



A(r) = A(r ) 



(8) 



where S is the expansion index, which is super-radial (S > 2) 
at the wind basis up to the merging radius when S becomes 
2. S is determined from Equation (7). 



2.2 The energy equation 

Falceta-Gongalves & Jatenco-Pereira (2002) showed that, 
even though presenting low effective temperatures, the tem- 
perature gradients at the wind basis of late-type supergiant 
stars could play an important role in accelerating the gas, 
as occurs in the Sun. In their model, the temperature profile 
was assumed to be a r dependent function, obtained from ob- 
servational data. The wind temperatures in giant cool stars 
are expected to increase from the photosferic value up to 
10 4 K, much lower than what is observed for the Sun, but 
still important in the wind acceleration if the gradient oc- 
curs in small lengthscales. From their calculations, Falceta- 
Gongalves & Jatenco-Pereira (2002) also showed that the 
thermal pressure is even more important than radiation 
pressure at the non-pulsating phase of these objects. 

In a consistent model, to avoid assuming any empirical 
function for the magnetic field geometry and to determine 
the wind temperature at each wind position (r), we have 
to solve the energy equation, which is determined from the 
balance between wave heating and the adiabatic expansion 
and radiative coolings (Hartmann, Edwards & Avrett 1982; 
Vidotto & Jatenco-Pereira 2006). Thus, neglecting conduc- 
tion, we write the energy equation as: 
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where e = ((SB) 2 ) / (4tt) is the wave energy density, which 



is described below (Sec. 2.3), Q is the wave heating rate, i.e. 
the rate at which the gas is being heated due to dissipation 
of wave energy, and Pr is the radiative cooling rate, both in 
erg cm~ 3 s _1 . The wave heating can be written as: 



Q = y(u + Va) 

and the radiative cooling is given by: 
Pr = A n e n H , 



(10) 



(11) 



where n e is the electron density, nu is the hydrogen den- 
sity and A is the radiative loss function. Here, we adopt 
the A function given by Schmutzler & Tscharnuter (1993) 
and calculate n e with the modified Saha equation given by 
Hartmann & MacGregor (1980). 



2.3 Wave energy density 

The wave energy density at each step may be calculated, 
using a WKB aproximation, from the wave action conser- 
vation. Under this assumption, the wave energy density is 
dissipated as follows: 



e = eo 
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where M — u/va is the Alfven-Mach number, va = 
(B/^/4np) the Alfven speed and L the wave damping length. 
Also, the wave flux (4>a) at ro is evaluated by (Jatenco- 
Pereira & Opher 1989): 



6 a = e Q v A o (l + ^Mo^J ■ 



(13) 



Actually, the WKB aproximation is not valid when 
the perturbation wavelength is much larger than the gra- 
dients lengthscales of the system, e.g. when A > L. As 
a consequence, part of the wave flux is reflected and the 
Alfven waves propagating on opposite directions decay as 
the generated beat waves interact with the gas particles. 
Davila (1985) pointed out that, in the Sun, the gradients 
of the coronal parameters are high enough to make the 
WKB aproximation not applicable. However, Usmanov et 
al. (2000) found good agreement between the observational 
data and models with this approach. In this model, we used 
Equation (12) assuming then linear perturbations, however 
included non-linear effects in the damping length function 
(Lnl). The non-linear damping length is given by (Jatenco- 
Pereira & Opher 1989): 

L„l -L (^Y { P^(1 + M), (14) 

\VAoJ {(Sv) 2 ) 

where ((Sv) 2 } is the averaged squared perturbation velocity 
amplitude and Lo is the damping length at the wind basis, 
which is mainly dependent on the assumed wave frequency 
spectrum (Lagage & Cesarsky 1983). Here, we will let it as 
a free parameter. 



2.4 The wind equations 

Finally, once the magnetic field structure is determined, we 
evaluate the temperature and velocity of the wind along the 
magnetic field. Using the Equation (8) and considering the 
magnetic flux conservation, Equations (1), (2) and (9) in the 
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Simplified geometry 
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Figure 1. The magnetic field structure for a constant super- 
radial index 5 = 5 (dashed line) and that determined self- 
consistcntly in our model (solid line). 
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Figure 2. Wind velocity profiles obtained for the initial damping 
lengths L = 0.2 ro (solid line), L = 1.0 ro (dotted line) and 
L = 5.0 ro (dashed line) as function of the distance. 



radial direction are reduced to (Vidotto & Jatenco-Pereira 
2006): 
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Equations (13) — (16) fully describe the wind param- 
eters and the magnetic field geometry under the given as- 
sumptions. In the next section, we show the main results by 
applying these equations in a typical cool supergiant star 
and compare them with previous works. 



3 RESULTS AND DISCUSSIONS 

For the last decades, the validity of a wind model was con- 
strained in reproducing only both the terminal velocity and 
the mass-loss rate of a given star. As a consequence, a num- 
ber of accelerating mechanisms were found in accordance 
to observations of these constrains. However, with the high 
resolution observations and more sensitive intruments the 
parameters radial profiles will become measurable and will 
be decisive on the modelling choice. Here we determine the 
velocity, density and temperature profiles and discuss their 
dependence on initial assumptions. 

We applied the described model on a cool supergiant 
star with M* = 16 M , r = 400 R , p = 10" 13 g cm -3 , 
Bo = 10 G and T = 3500 K. We also assumed a filling 
factor a = 0.1, according to solar observations. 



3.1 Wind profiles 

To study the importance of the damping length on the wind 
profiles we performed calculations setting the initial wave 



flux as (/>ao — 3.10 7 erg cm~ 2 s _1 , and varied Lo as 0.2, 1.0 
and 5.0 ro. The magnetic field structure is shown in Figure 
(1). The solid line represents the solution of Equation (7) for 
the given stellar parameters, compared to a simple constant 
expansion factor 5 = 5 (dashed line), as used by Falceta- 
Gongalves & Jatenco-Pereira (2002) . Comparing both curves 
we find that the self-consistent geometry is more divergent 
near stellar surface, and the flux tubes are expected to merge 
at lower distances, in agreement with more refined models 
for the Sun (Esser et. al. 2005). We obtained an initial value 
of S(ro) ~ 35, that decreases with distance to the stellar 
surface. From the assumed a, we found that the neighbour 
fluxtubes merge at r = 1.41 ro. At this point, with a sudden 
decrease of S from its super-radial value to S = 2 (radial), 
the wind properties change and it is noticeable as breaks in 
the temperature, density and velocity profiles. 

We obtained the same geometry for the following calcu- 
lations considering different initial parameters. This because 
the super-radial magnetic field geometry occurs mainly at 
the wind basis, where the magnetic pressure is highly domi- 
nant. In this sense, we expect the expansion index (S) to be 
independent on the other initial parameters such as the wind 
velocity, the wave flux densities and the damping length. 

In Figure (2) we compare the wind velocity profiles for 
the different initial damping lengths. As expected, the larger 
the damping length the higher is the wind velocity. However, 
for previous works that do not take into account any high 
initial divergence, the wind velocity is even larger. The ob- 
tained divergence result in a fast dilution of the wave energy 
density near the stellar surface. As a consequence, the lack 
of wave flux at larger distances results in lower velocities. 
This result shows that even low damped wave fluxes can 
be responsible for low terminal velocities, depending on the 
field divergence. 

In Figures (3) — (5) we present the gas mass density, 
wave energy density and the temperature, respectively, for 
the same calculations. For the gas density, all three lines 
show the same behaviour at r < 1.5 ro indicating that the 
flux tube divergence is the dominant factor if comparing to 
the gas expansion due to acceleration. This is an interest- 
ing result since, typically, the density profile is believed to 
be closely related to the velocity profile only. In Figure (4), 
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Figure 3. Wind density profiles obtained for the initial damping 
lengths L = 0.2 ro (solid line), L = 1.0 ro (dotted line) and 
L = 5.0 ro (dashed line) as function of the distance. 
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Figure 5. Wind temperature profiles obtained for the initial 
damping lengths L = 0.2 ro (solid line), L = 1.0 ro (dotted line) 
and L = 5.0 ro (dashed line) as function of the distance. 
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Figure 4. Wave energy density profiles obtained for the initial 
damping lengths L = 0.2 ro (solid line), L = 1.0 ro (dotted line) 
and L = 5.0 ro (dashed line) as function of the distance. 

we notice that, in despite of what was obtained in previ- 
ous works, the wave energy density is weakly dependent on 
the initial damping length at the wind basis. This mainly 
because, again, the field divergence is dominant on the di- 
lution process. On the other hand, for distances larger than 
the flux tubes merging position the wave damping is domi- 
nant. In Figure (5), the temperature profiles show different 
features. For the given stellar parameters, if Lo = 5.0 ro, the 
wave heating is less effective and the radiative losses added 
to the expansion cooling result in a negative temperature 
gradient, which occurs in a very narrow region until the 
density becomes low enough for the radiative losses become 
less important. For highly damped waves, the wave heating 
dominates the radiative losses and the temperature gradient 
is positive near the surface. Another interesting feature is 
the maximum temperature position. Since the temperature 
gradient is mainly dependent on the wave damping (via Q 
parameter), we notice that the lower is the damping length, 
the closer to the stellar surface the maximum temperature 
occurs. 

Recently, Zeeman splitting of H2O masers observations 
led Vlemmings, Diamong & van Langevelde (2002) to in- 
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Figure 6. Wind velocity profiles obtained for initial magnetic 
field strengths B = 10 G (dotted line) and B = 100 G (solid 
line) as function of the distance. 



fer magnetic field strengths on supergiant stellar surfaces of 
Bo > 100 G. Also, theoretical dynamo models for AGB stars 
reveal that the planetary nebulae shapes can be explained 
by magnetic field strengths of ~ 200 G (Blackman et. al. 
2001). For this reason, we also perfomed numerical calcula- 
tions for Bo = 10 and 100 G, setting the initial parameters 
for a low damped wave flux (Lo = 5 ro). We also assumed 
the same wave amplitude relative to the magnetic field for 
both cases. 

The velocity profile for both initial magnetic field 
strengths are shown in Figure (6). For the larger magnetic 
field strength (Bo = 100 G) we obtained the higher termi- 
nal velocity (u > 200 km s -1 ). This effect is mainly due to 
the higher wave energy flux assumed at the wind basis to 
accomplish the same relative wave amplitude. 

The temperature profiles, shown in Figure (7), also 
present features dependent on the magnetic field strength. 
As for the velocity, a high wave energy flux results in a 
high temperature, at least near the stellar surface due to 
the higher energy density. 
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Figure 7. Wind temperature profiles obtained for initial mag- 
netic field strengths Bo = 10 G (dotted line) and Bo = 100 G 
(solid line) as function of the distance. 

3.2 Best-fitting model 

For stars with parameters similar to those used in the pre- 
vious subsection, observational data reveal typical mass-loss 
rates of M ~ 1CP 7 — 1CP 6 Mq yr _1 and terminal velocities 
of u ~ 70 km s _1 . Unfortunately, the available data are lim- 
ited in spatial resolution for most of the stars, and it is not 
possible to fit the complete radial profiles. 

Assuming a surface magnetic field strength B = 10 G, 
and a low damped wave flux (L = 5 ro), it was possible 
to reproduce both the wind terminal velocity and the mass 
loss rate using an Alfven waves flux of 4>ao = W 7 erg cm~ 2 
s _1 at the wind basis. This value corresponds to a wave 
amplitude of {(SB) 2 }? ~ 3.10~ 2 B , which is very plausible 
for a turbulent medium as that at the stellar surface. 

The velocity and temperature profiles for this case are 
shown in Figures (8) and (9), respectively. The velocity pro- 
file reveals a peak of u > 100 km s _1 at r < 2.0 ro, and 
slightly decreases for larger distances until reaching the ob- 
served value. The temperature profile presents an initial neg- 
ative gradient reaching temperatures T < 2500 K in a nar- 
row region. Also, near r = 1.5 ro the temperature reaches 
the maximum value of ~ 6000 K. For r > 3.0 ro, where the 
wave heating and the radiative losses are low, the tempera- 
ture decreases due to the adiabatic expansion. 



4 CONCLUSIONS 

We propose a self-consistent wind model to determine the 
parameters profiles for a supergiant late-type star. To de- 
termine the magnetic field geometry we used an expansion 
method over the wind physical parameters as proposed by 
Pneuman, Solanki & Stenflo (1986). Near surface, the mag- 
netic pressure inside the flux tubes are higher than that of 
the surrounding medium, forcing the field lines to curve. We 
found an initial super-radial expansion factor S > 30 at the 
wind basis, much higher than the value used by previous 
authors that included empirical relations to account for the 
magnetic field geometry based on solar observations. 

Considering a supergiant late-type star with M* = 16 
M Q , r = 400 R , p = 10~ 13 g cm~ 3 , B = 10 G and 




r(r ) 

Figure 8. The wind velocity profile for the best fitting parame- 
ters in the case of low damped waves. Here we assumed Bo = 10 
G, Lq = 5.0 ro and <f>A0 = 10 7 erg cm~ 2 s _1 . 
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Figure 9. The wind temperature profile for the best fitting pa- 
rameters in the case of low damped waves. Here we assumed 
B = 10 G, L = 5.0 r and <pA0 = 10 7 erg cm -2 s _1 . 

To = 3500 K, we obtained the wind velocity, density and 
temperature profiles. Typically, Alfven wave driven winds 
result in high velocity winds (u > 100 km s _1 ), unless some 
strong wave damping mechanism takes place at the wind 
basis. We showed that this conclusion is correct in the case 
of low divergent magnetic field structures. In this work, the 
strong divergence is responsible for a rapid wave spatial di- 
lution near surface, resulting in a lower wind velocity even 
for low damped waves. 

We reproduced both the typical mass-loss rate (M ~ 
10~ 7 - 10~ 6 M Q yr _1 ) and terminal velocity (it ~ 70 km 
s _1 ) observed for these objects, assuming a weakly damped 
(Lo = 5.0 ro) Alfven wave flux of 4>ao = 10 7 erg cm~ 2 
s _1 . The velocity profile reveals an efficient acceleration at 
r < 1.5 ro, reaching the maximum value ~ 100 km s _1 . In 
this region the wind is mainly accelerated by the wave energy 
density and the thermal pressure gradients. Afterwards, the 
absence of the wave acceleration and the cooling gas result 
in a decrease of the velocity to the observed values. For 
the temperature, assuming a weakly damped wave flux, the 
radiative losses and the expansion cooling are dominant near 
surface, and the temperature gradient is initially negative. 
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The temperature falls to ~ 2500 K in a sharp region and 
then, as density decreases as the wind accelerates and the 
flux tube expands, it increases up to ~ 6000 K at r < 2.0 
ro. For higher distances, where the radiative losses are low 
and the wave heating is no longer effective, the temperature 
decreases mainly due to the adiabatic expansion. 
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